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Abstract 



Recently, various cross sections of e + e~ annihilation into hadrons were accurately 
measured in the energy range from 0.37 to 1.39 GeV with the CMD-2 detector at the 
VEPP-2M collider. In the 7r + 7r~ channel a systematic uncertainty of 0.6% has been 
achieved. A Monte-Carlo Generator Photon Jets (MCGPJ) was developed to simulate 
events of the Bhabha scattering as well as production of two charged pions, kaons and 
muons. Based on the formalism of Structure Functions, the leading logarithmic contri- 
^ ' butions related to the emission of photon jets in the collinear region are incorporated 

into the MC generator. Radiative corrections (RC) in the first order of a are accounted 
Qh! for exactly. The theoretical precision of the cross sections with RC is estimated to be 

better than 0.2%. Numerous tests of the program as well as comparison with other MC 
generators and CMD-2 experimental data are presented. 



x. 

c3 . 1 Introduction 

The cross sections of e + e~ annihilation into hadrons are very important in various problems 
of particle physics and, in particular, they are required for the evaluation of the hadronic 
contribution to the anomalous magnetic moment of muon, = (g — 2)^/2. The recent 
measurement of a M at BNL pQ led to a new world average differing by 2.7 standard deviations 
from its theoretical evaluation. One of the main ingredients in the theoretical prediction 
for a M is the hadronic contribution related via a dispersion integral to the cross section of 
e + e~ annihilation into hadrons. In the case of a^ ad , the VEPP-2M energy range gives the 
major contribution both to the hadronic vacuum polarization contribution itself and to its 
uncertainty [3 E] . This uncertainty is dominated by systematic errors of the experimental 
values of R(s) which are used as an input to the integral with the proper kernel function [I]: 
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The quantity R(s) is denned as R(s) = a(e + e~ — > hadrons)/a(e + e~ — > /i + and at high 
energies can be calculated within the QCD framework whereas at low energies the experimental 
data are required. A numerical computation of this integral can be found elsewhere j2] and in 
relative unities its evaluation yields the result ~ 70 ppm. 

The goal of the new BNL experiment [5] is to measure the anomalous magnetic moment 
of muon with the relative accuracy ~ 0.25 ppm. To reduce the current systematic error of the 
hadronic contribution to a^ ad at least to the same level, the theoretical precision of the cross 
sections with radiative corrections (RC) should be better than 0.3% as it follows from a simple 
estimation: 70 ppm xO.3% ~ 0.2 ppm. This short observation shows why the knowledge of 
the cross sections e + e~ annihilation into hadrons with high precision is extremely important. 

The detection efficiency, background conditions and criteria of event selection from the raw 
data differ for specific e + e~ annihilation modes. Therefore, expressions for the cross sections 
with RC on which the M C generator is based, should have a completely differential form with 
respect to the kinematic variables of final particles. In this case the influence of the selection 
criteria as well as the trigger efficiency and many other specific resolutions of the detector can 
be naturally incorporated in a MC generator. 

During the last thirty years considerable efforts were devoted to elucidate theoretical un- 
derstanding of the accuracy of cross sections with RC, particularly in the case of e + e~ and 
7T + 7r~ pair production at low energies. The radiatively corrected cross sections for annihilation 
channels with an accuracy of about 0.1% were obtained in 6\. Unfortunately, expressions for 
these cross sections do not contain the angular distributions for the emitted photons and, as 
a result, it is not possible to reconstruct the kinematics of the final particles correctly. On the 
other hand, the differential cross sections were obtained in [Jj, but their relative accuracy is 
about 1%, since only 0(a) QED corrections were taken into account. 

The work [S] is based in part on a combination of the approaches [SJ Ej mentioned above. 
To achieve the accuracy ~ 0.2%, higher order radiative corrections were taken into account by 
means of the Structure Function (SF) formalism jH]. It involves a convolution of the boosted 
Born cross section with the electron (positron) SF, which describes the leading effects due to 
emission of photons in the collinear region as well as radiation of e + e~ pairs. These enhanced 
contributions are proportional to (a/tr) 11 ln n (s/m 2 ,), n = 1,2,... and are referred to as the 
leading ones. Moreover, in the smoothed representations of the SF [S| a certain part of these 
corrections is exponentiated and evaluated in all powers of n. The non-leading contributions 
proportional to (ex /it) are incorporated exactly according to [6 3 by means of a so-called K- 
factor. The next-to-leading contributions of the second order (a / n) 2 \n(s / m 2 ,) ~ 0.01% are 
fortunately small and can be omitted, keeping in mind the intended precision tag 0.2%. 

The vacuum polarization effects in the photon propagator are treated as in [H] for the 
lepton channels. These effects are not included in RC for the hadronic modes according to 
the generally accepted agreement ^Uj- The emission of one hard photon at a large angle is 
described by a differential formula, which allows to take into account specific experimental 
conditions and cuts. Based on numerical calculations, which will be given below, it was 
established that the O(o) radiative corrections together with photon jet radiation in a collinear 
region are sufficient to achieve a relative precision of the cross sections about ~ 0.2%. 

The purpose of this paper is to describe the Monte-Carlo Generator Photon Jets (MCGP J) 
which simulates processes e + e~ — > e + e~, 7r + 7r~, K + K~ and K^Ks- This generator was 

used while CMD-2 experimental data were processed. The MCGPJ code has a modular struc- 
ture that simplifies the implementation of additional hadronic modes as well as the replacement 
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of matrix elements of the current cross sections by a new one. The effects of the final state ra- 
diation (FSR) for the channels tt + tt~ , K + K~ have been incorporated into the program. 
The pions were assumed to be point-like objects, and the scalar QED was applied to calculate 
virtual, soft and hard photon emission by charged pions (kaons). 

The relevant formulae for the cross sections with RC of order a are collected here from 
many other papers. It is done specially to have a possibility of quantifying the difference 
between cross sections due to radiation of one photon and photon jets in the collinear region. 
On the other hand, some expressions for RC are revisited in the present paper and explicit 
analytical formulae will be presented in the form convenient for the MC generator construction. 

2 Monte-Carlo generator for events of Bhabha 
scattering at large angles 

The boosted Born cross section of the process e~ + e + (z 2 p+) — ► e ~ (pi) + e + (P2) , corrected 

for vacuum polarization factors in the s and t channels, when initial particles lose some energy 
by radiation of photon jets in the collinear region, has the following form jH] in the cm. frame: 

d ~e+e-^ e +e- ^ ^ ^ 4^ ^2 ^ ~2 + ~2 p + ^2 



dfii a 2 s \2t 2 \l -II(t)| 2 2§ 2 |l-n(s)| 2 

+ MiF(i-n( s ~))(i-n(f))})' (1) 

where z\ and z 2 are the electron and positron fraction energies after radiation of photon jets 
(-21,2 = £1,2/ ' £ beam) j II (s) and n(t) are the vacuum polarization operators in photon propagators 
in the s and t channels, respectively. The Mandelstam variables in the Lab and c.m.s. are 
defined as usual: s = 2p_p + , t = — 2p-pi, u = —2p_p 2 , s = sz\z 2 , t = —sz\Yi(l — ci)/2, 
u = -sz 2 Y 1 (l + ci)/2, si = 2p 1 p 2 , ti = -2p + p 2 , ui = -2p + p 2 , c\ = cos8 1 , where 9i is a 
polar angle of the final electron with respect to the electron beam direction, Y\ and Y 2 are 
the relative energies of final e~ and e + . The energy-momentum conservation law allows to 
reconstruct the kinematics of final particles and to find these energies and a positron polar 
angle 6 2 : z\ + z 2 = Y\ + Y 2 - energy conservation; z 1 — z 2 = Y\ cos6\ + Y 2 cos9 2 - momentum 
conservation along the Z-axis; Y\ sin d\ = Y 2 sin 9 2 - momentum conservation in the plane 
perpendicular to the Z-axis. From these equations one can find that 

2 Zl z 2 (zf + zp- (z 2 - z\)c x 

a a 

(zf — z\) — (z\ + z 2 )ci , 
c 2 = 2 2 7-2 -2T— , where a = z 1 + z 2 -(z 1 - z 2 ) Cl . (2) 

\ Z 1 + Z 2) - \ Z l - Z 2) C l 

The expression for the differential cross section with one photon emission in the reaction 
e (P-) + e+ (p+) ~^ e (Pi) + e+ (P2) + j{k), was obtained in J\ (see also references therein) 
and reads 

j e+e-^e+e"7 °^ r>e+ e~ ^e+e~ ^Pl ^P2 

d ^hard = 7T^- R hard Z ^ [P- + P+ ~ Pi ~ P2 ~ ty, (3) 

where ei, e 2 , and uj are the energies of the final state electron, positron and photon, respec- 
tively; 5-function provides the energy-momentum conservation. 
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The expression for the quantity R\ard ^ e e 7 which contains the vacuum polarization effects 
in photon propagators was derived in [HJ: 



pe 1 e 
^hard 



+e T e 7 



(WT) n 
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m 2 /s 2 + (s + £) 2 t 2 + {s + t) 2 f (s + tf 



+ 3?e 



x' + 2 V2t 2 (f-n(t)) 2 2s 2 |i-n( s )| 2 U(i-n(s))(i-n(t)) 

m 2 e f s 2 + (s + t 1 ) 2 | tj + js + h) 2 | (s + t^ 2 



x'_ 2 V2t 2 (i-n(t 1 )) 2 2s 2 |i - n( s )| 2 Ui(i-n(s))(i-n(*i)) 

m 2 e ( sj + (si + t) 2 t 2 + (si+t) 2 f (Si+t) 2 



i U 2 (i - n(t))2 + 2s 2 |i - n( Sl )l 2 W(i - n( Sl ))(i - n(t)) 



"x 2 V2t 2 (i-n(t!)) 2 1 2s 2 |i-n( Sl )| 2 1 ^U 1 t 1 (i-n( Sl ))(i-n(t 1 )) 

where x± — kp± and x± = &Pi,2- The quantity (WT)n describes the process with one hard 
photon emission and gives the dominant contribution outside the collinear region [HJ: 

(WT)u = I -. J! , , + SlSl TT 



|i - Ti(s)\isx'-x' + |i - n( Sl )l 2 ^x-x+ |i - iWtx+xV 



TlTl + Ke 



li-n^oi^x-x'- 



(i-n(t))(i-n(* 1 ))tt lX -x , -x+x+ 



551 TS* 

+ 



(f - n(s))(i - u(s 1 ))*ss l x-x-x+x+ (i - n(t))(i - n( s ))ts X '-x+x' + 

T\ S i T± S 



(i - n(tO)(i - n( Sl ))t lSlX _ x '_ x+ (f - n(t x ))(i - n(s))t 1 s X ^x'-x' + 

TSi 



(i-n(t))(i-u(s 1 ))ts lX -x+x' + 

where the following notations were used: 



(5) 



SS = SxSt = t 2 + tj + u 2 + uj, 
TT = TxTx = s 2 + sj + u 2 + uj, 

SSi = if + tj + u 2 + uj) x (tx+X+ + ^lX-X- - «X+X- - «iX-X+), 
TTi = (s 2 + s 2 + u 2 + uj) x (mx+X- + «iX-X+ + SX-X+ + siX-X+), 

T5 = -^(M 2 + « 2 )(s(t + Si) + t(s + *i) - UUi), 

TSi = -^(u 2 + uj){t( Sl + ti) + si(s + t) - uui), 
TiS = ^{u 2 + w 2 )(ti(s + t) + s(si + h) - UUi), 

TiSi = l -{u 2 + uj)(si(s + tO + ti(si +t)- uui). (6) 

The main contribution to the cross section due to photon radiation comes from the collinear 
region, where the cross section exhibits very steep behavior JT] . The collinear region is a part 
of the angular phase-space with four narrow cones (Fig.[TJ surrounding the directions of motion 
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of the initial and final particles. The emitted photon should be inside these cones with an 
angle 29 . The angle 9 should obey to the follow restrictions, I/7 <C 9 <C 1, where 7 = e/m e . 
It serves as an auxiliary parameter and usually its value is taken at about ~ 1/ ^fty. The cross 
section integrated inside these narrow cones according to [S] is: 



d(T coll 



+ 



dcr] 



e 1 e 




dOi 





2i 



X 2 \ / 1 , 0g\ X 

z H L-l+k-rM 

2 4/2 



(7) 



where L = ]n(s/m%), z = 1 — x and the boosted Born cross section is defined in Eq. (pQ). The 
auxiliary parameter A =Ae/e (A <C 1) serves as a separator of hard and soft photons, e is 
the beam energy. The terms proportional to (a/7c)(L — 1) are accounted for in the SF jH] and 
therefore should be removed from this expression to eliminate the double counting. 

The remaining four terms can be interpreted as the so-called compensators which cancel 
out the dependence of the total cross section on the auxiliary parameter # when they are 
summed with the cross section Eq.fjHJ) describing one hard photon emission outside cones. 

Collecting all the discussed above terms into one formula we get the complete expression 
for the master formula describing the process e + e~ —>■ e + e~ + nj, which can be presented as: 



^0-e + e — >e + e +ri~f 
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dxij ix 2 J ix 3 J dxd 
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21 daf e " 
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(23 + =2) In ^ + 
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A 
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d#j 



7T 
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In - In A H — 

t 2vr 2 s 



dfii 



1,1: 



6 (cuts) 



fe°>A £ 



wr e+e " 7 ^0(cuts), 



(8) 



where £1,2,3,4 are the relative energies of photon jets emitted along the initial and final particles; 
21,2,3,4 = 1 — £1,2,3,4 are the energy fractions of electrons and positrons after radiation of photon 
jets; 6 (cuts) is a step- function equal to 1 (0) if the kinematics variables obey (or not) the 
selection criteria; the expression for Ksvifii) can be found in 0|H]. More details concerning 
the implementation of the SF formalism as adopted in the present paper can be found in j^j. 
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The SF approach provides essential improvement of accuracy for the Bhabha cross section 
by taking into account radiation of photon jets in the collinear region. These improvements 
as well as others performed in [8] are summarized below: 

1. The radiation of photon jets (enhanced contributions) is taken into account by means of 
the SF formalism. 

2. To combine the cross sections, describing radiation of one hard photon inside and outside 
narrow cones, the four compensators are embedded into the master formula (EqIBJ) • 

3. The boosted Born cross section contributes to the total cross section in conformity with SF 
weights in the master formula (EqJSJ) . 

4. The vacuum polarization effects inserted into all photon propagators exactly. 

5. Non-leading contributions of order a are accounted for by means of the so-called K-f actor. 
The integration limits of the first term in EqlHlwere divided in two parts from to Ae and 

from Ae to the maximal jet energy. As a result, the four-fold integral splits into sixteen parts. 
Those of them with one photon jet radiation are merged in a proper way with compensators 
in the master formula. In this case, the total cross section is subdivided into seventeen cross 
sections with own specific kinematics - number of photon jets and directions of their radiation. 
The first contribution accounting for effects due to soft and virtual radiative corrections is given 
by 



da{ 



+ e — >e+e +n-y 



A A A A 



1 1 J J dxidx2dx 3 dx4V(zi,s)V(z 2 ,s)V(z3,s)T>(z4,s) 





x ( i + °^) d95 T'^'-^^ihA^ (M) . 

7T aili H \t J dili 

The photon jet energy emitted by each charged particle can be up to Ae. This part also 
contains the contribution due to production of virtual and soft real e + e~ pairs if 2m e < Ae. 

The next four terms represent the contribution to the cross section with one hard jet 
emission along the motion of any charged particle, supplied with the virtual and soft leading 
logarithmic corrections. One of these terms with the relevant compensator is: 

d(j e +e -_> e+e - +n7 1AAA ,__ + 



III J dxidx 2 dx 3 dx4V(z 2 , s)V(z 3 , s)V(z 4 , s) — — 



dQi J J J J ™ ° * ' v ' ' dQi 

A 



a ~, s a 1 / , xi , , 0n x 



V{z x , s)(l + -K SV ) + -- l( Zl + ^) In J + ^ 

7T 7T Xi \ 2 4 2 



e(cuts). (io) 



The other similar terms can be written out in the same way by the permutation of limits 
between integrals. 

The next six terms represent the contribution to the cross section when two jets are emitted 
simultaneously along momenta of any two charged particles. One of these terms reads 



dVli 



dxidx2dx 3 dx i V{zi 1 s)V(z 2 , s)V(z 3 , s)V(z / i, s) 



a a o o 



dtti \ 7T 



1 + -K SV e(cuts). (11) 
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The other similar terms have the identical structure and are obtained by the limits permuta- 
tion. 

The following four terms represent the contribution to the cross section when three photon 
jets are emitted along momenta of any three charged particles. One of these terms is given by 

j e + e — — >e+e _ +n7 \ \ \ Q 



III f dxidx 2 dx3dx4 : V(z 1 ,s)V(z2,s)V(z3,s)V(z4,s) 



A A A 



dag e ^ e e (z 1 ,z 2 ) , a r> 



x v " z/ (^1 + -K sv j e(cuts). (12) 

The cross section with emission of four jets along the momenta of each initial and final 
particle is written below, 



e + e — >e+e +wy ' ' ' 



III j dxidx2dx 3 dx4V(zi, s )'^ > ( z 2, s)V(z 3 , s)V(zi, s) 



A A A A 



dog— 2 )r 1 + ^ 5v) e (cu ts). (is) 



df^x \ 7T 

The cross section with one hard photon emission outside the collinear region reads 
dv e i7 e ^ e+e 7 « 3 f De + e -^ e + e -+ 7 SXiXdxdQ^ 

An = / R hard W\ ■ 2 / /0 N e ( cuts )> 14 

di 2i ivr^s 7 8(1 — ism w 2) 

s 7 >e 

where ^ is an angle between the momenta directions of photon and final electron. 

The cross section calculation was performed by the Monte-Carlo method. The cutoff 
energy Ae was chosen at ten electron masses to optimise the efficiency of event simulation 
(A = Ae/e ~ 1%). Since the master formula is singular on some variables, the main of them 
have been isolated: photon energy and emission angle were generated according to functions 
\/u)(e — u) and 1/(1 — /3| cos 2 # 7 ), respectively. The main contribution to the Bhabha cross 
section comes from the t channel and it was generated by the function 1/(1 — cos#i) 2 . 

The selection criteria adopted here for the simulated events are similar to those used in 
CMD-2 data analysis jHj and they are: 

the acollinearity cut in the scattering plane is \A8\ < 0.25 rad, where A# = 9\ + 62 — 71; 
the same for azimuthal plane, \A(f)\ < 0.15 rad, where A<p = \<pi — 02 1 — 7r; 
the angular acceptance is 1.1 < avCT < n — 1.1, where aver = {6\ — 9 2 + 7r)/2; p± 2 > 90 MeV/c. 
Below, if nothing specially told about selection criteria the latter will be taken in mind. 

The MCGPJ program consists of two main stages. In the first stage with a soft selection 
criteria all majorants for seventeen parts are determined, in the second one the cross sections 
with experimental selection criteria are being determined. The generator simulates an event 
according to specific kinematics for each cross section. The weight of event is determined as the 
ratio of a given cross section (one from seventeen) to the total one. The kinematic parameters 
of the simulated events are stored in the proper histograms, which can be compared with 
experimental distributions. 
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Numerous tests have been performed for the c.m.s. energy of 900 MeV. The cross section 
dependence on the auxiliary parameter Ae is shown in Fig|2] after integration over the re- 
maining kinematic variables. It is seen that the cross section variations are inside the claimed 
precision while Ae changes by a factor of 10 4 . The cross section variations with an auxiliary 
parameter 8 do not exceed ±0.1% level (FigJ2J). From that certainly follows that the cross 
section stability to the auxiliary parameters Ae and at least is not worse than 0.1%. 

Comparison of different kinematic distributions simulated by the MCGPJ generator and 
BHWIDE ^2] was performed. Distributions over parameters AO = Q\ + 82 — it and A</> = 
101 — <^>2 1 — 7T are plotted in FigslH El Good agreement is seen while AO and A0 vary in the 
wide range. Asymmetric shape of the distribution on A0 is due to specific kinematic of the 
events, when photon and e + e~ pair fly in diametrically opposed sides. The cross section for 
such kinematic has enhancement which visualizes as shape asymmetry of the distribution. 

The event distributions produced by both generators are presented in Fig. El as a function 
of missing energy (e m i S , = 2e — S\ — e<i\ Both distributions are close to each other except for 
the energy region where soft and hard photons are merged. A visible bump is observed in this 
point. This bump originates from a slightly different dependence on the cutoff energy of both 
the compensators and the cross section with one hard photon. The cut on acollinearity AO ~ 
0.25 rad, is equivalent to missing energy ~ 100 MeV that is rather far from cutoff energy Ae ~ 
7 MeV. As a result, the contribution into the total cross section connected with this spurious 
features is negligible. 

The relative difference of the cross sections, presented in FigJTJ calculated by the MCGPJ 
code and BHWIDE is less than 0.1% for the VEPP-2M energy range. This difference versus 
the acollinearity angle AO is plotted in FiglHl One can see that the size and sign of the 
difference depend on the particular choice of AO. The reason of the difference about ~ 0.5% 
for AO ~ 0.05 rad arises from the fact that all photons (except one in our code) are emitted 
strongly along the motion of electrons (positrons) whereas in BHWIDE they have some angular 
distribution. The difference of ~ 0.3% for the large acollinearity angles \A0\ ~ 1 rad is due 
to the fact that the BHWIDE code simulates one hard photon only. It is worth noticing the 
MCGPJ code describes the shape of tails of the different kinematic distributions a bit correctly 
and, as a result, it is preferable for applications when soft selection criteria are used. 

It is important to reliably estimate the total theoretical precision of this approach. In 
order to quantify a theoretical error, the independent comparison has been performed with 
the generator based on Ref. where 0(a) QED corrections are treated exactly. It was found 
that the relative difference of cross sections is more than 1% for small acollinearity angles 
AO < 0.1 rad ( Fig. El ) and it is less than ~ 0.2% for acollinearity angles ~ 0.25 rad. From 
that immediately follows that the radiation of two and more photons (jets) in the collinear 
region contributes to the cross section by the amount ~ 0.2% only. Therefore, we can conclude 
that the theoretical precision of the Bhabha cross section with RC is certainly better than 
~ 0.2% for our selection criteria. 

The EM calorimeter of the CMD-2 detector allows to separate Bhabha scattering events 
from others with high confidence level The distributions in acollinearity angles AO and 
Acf) are presented in Figs llOlfTl To increase the experimental statistics, all the CMD-2 data 
at energies greater than 1040 MeV collected in these plots. The momentum and angular 
resolutions, interaction with the detector material were added to the kinematic parameters 
of simulated events. The histograms were fitted by two Gaussian functions. Their relative 
weights and widths were fit parameters. Good agreement between experiment and simulation 
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can be seen. 

The agreement between experiment and simulation becomes significantly worse when the 
MC generator based on Ref. [7] with 0(a) corrections is used. It is seen in Figs ll2lf3l where 
two-dimensional plots are presented. The points in these plots correspond to the electron and 
positron energies. Different population of events is observed far aside from the area where 
semi-elastic events are concentrated. About ~ 1% events have correlated low energies and 
they are distributed predominantly along a corridor which extends from the right upper to 
the left bottom corner of this plot. The appearance of these events is due to simultaneous 
radiation of two jets with close energies along either initial or final particles. The condition 
pj; 2 >90 MeV/c is very soft and only owing to this fact the integrated cross sections are equal 
to each other within ~ 0.2%. If p± 2 is about ~ 220 MeV/c, the relative difference increases up 
to ~ 1% (FigJ14jl. For the values p\ 2 more than 350 MeV/c the difference changes a sign and 
grows up. The cross section with photon jets becomes smaller than with one photon under 
condition p± 2 > 350 MeV/c. This feature has a simple explanation. The distribution width 
of semi- elastic events in the first plot is broader than for the second one due to radiation of 
many soft photons and, as a result, these events are smeared more broadly near the peak area. 



3 Monte-Carlo generator for production of muon pairs 

The same approach was used to create the MC generator to simulate production of muon pairs 
in the reaction e~(z\pJ) + e + (z 2 p + ) — ► ^ (pi) + /^(j^), when initial particles radiate some 
energy by emission of photon jets in the collinear region. According to jH] the boosted Born 
cross section modified by the vacuum polarization effects in the photon propagator reads 



daf e ~^ + »~( Zl ,z 2 ) _ a 2 1 yi[z 2 1 (Y 1 -y 1 c 1 ) 2 + z%(Y 1 + y 1 c 1 ) 2 + 8z 1 z 2 ml/s] 

dfii 4s | 1 - n(^iz 2 s) | 2 z\zl[zi + z 2 - Oi - z 2 )c 1 Y 1 /y 1 ] 



,(15) 



where y\ 2 — Y 2 2 — 4m 2 / s; Y"i j2 = Si, 2 /e are the muon relative energies; Zi }2 = 1 — xi,2, 
%i,2 — ^1,2/ £ are the relative energies of photon jets; c\ = cos#i, 9\ is a polar angle of 
negative muon. The energy-momentum conservation law, 



z 1 + z 2 = Y 1 + Y 2 , zi - z 2 = y x c x + y 2 c 2 , y\\j\- c{ = y 2 ^l - c£, 
allows to determine Yi, Y 2 and a positron polar angle 6 2 (c 2 = cos #2): 

1_ s z 1 z 2 + [z 2 1 z 2 2 -(ml/s)((z 1 + z 2 ) 2 ~(z 1 -z 2 ) 2 c 2 1 )]y 2 ^ z 1 + z 2 -c 1 (z 1 -z 2 y 1 ] 

The charge-even part of the cross section in the first order in a comes from one-loop 
virtual (V) and soft (S) radiative corrections and according to Ref. [TH] is given by: 

daf+r daf^' (1, 1) 2a + 



df2i dOi 7T 
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The expression for the quantity K£ ven was derived in 

'1 + /3 2 1 1 \ , 1+P 



even 



-1+P 

1-/3 2 



2(3 



1 1 , , 
- H + In • 

2 4/3, 



and reads 

( 1 1 + ^ 

V2/9 /? 



2/5 2 - /3 2 (1 - cf) 

. 1 + /3 .a 

In ^ s , p = In — 

1 — p m 1 



+ 



1 + (3 2 \k 2 
2/5 

L 



+ 2Li 2 



1+/3 



In 



+ 1/3 In 



r df 

U 2 (x) = -J T ln(l-0- (1? 



1 + /? 



2/? 2 



In the ultra relativistic limit ((3 — ► 1) the quantity takes the same form as for electrons 
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(19) 



The charge-odd part of the cross section comes from the interference of the Born amplitude 
with box-type diagrams as well as with amplitudes of soft photon emission by initial and final 
particles and is given by [HI |H] : 
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(20) 



For the ultra relativistic limit the same result as in 1-1] is obtained. 

The cross section of muon pair production with one hard photon emission is studied in detail 
elsewhere [T3*l |H] . This cross section in the differential form, keeping the relevant information 
about the kinematics of final particles, can be written as: 



a 
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(21) 
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where (3% is a velocity of negative muon. The quantity R^ d 7 consists of three terms 

and represents the cross section with one hard photon emitted by the initial and final particles 
as well as their interference: 
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Mandelstam variables and x±? X± are defined as for electrons. Similar to the Bhabha cross 
section, the master formula describing the process of muon pair production reads jH] 
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(23) 



where K = n 2 /6 — 1/4 + Kg Yen (s, 0\) + K^ dd (s, 61) and Q\ is a polar angle of negative muon in 
the cm. frame. The master formula drawn above provides within the scope of the discussion 
the intended cross section accuracy ~ 0.2%. The integration limits of the first term in Eql2*3l 
were divided into two parts as for electrons. A two-fold integral splits into four separate 
contributions. Those of them which describe one photon jet radiation are combined in a 
proper way with two compensators in the master formula. All other steps to construct the 
MC generator to simulate production of muon pair are similar to that as for electrons and can 
be found in [To] . 
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Numerical comparisons with the KKMC generator have been performed. The the- 
oretical accuracy of the formulae on which KKMC is based is about ~ 0.1%. The existing 
code in KKMC does not provide the correct description of vacuum polarization effects in the 
photon propagator at low energies, so they were switched off in both generators. The relative 
difference between cross sections calculated with the MCGPJ generator and KKMC in the 
VEPP-2M energy range is presented in FigIT5l Good agreement at the level of our precision 
±0.2% is observed. 

In the low energy range, the momentum resolution of the CMD-2 detector is sufficient to 
distinguish pions, muons and electrons. Thus we have a direct way to compare the number 
of selected muons to that of electrons divided by the ratio of the theoretical cross sections, 
a(ee — ► fifi)/a(ee — > ee) and check thereby the theoretical precision of the formulae with RC 
from experiment. The results for this double ratio are presented in Figllfil One can see that a 
deviation from unity of the double ratio does not exceed on average 1.4% with statistical and 
systematic errors about ~ 1.5% and ~ 0.7%, respectively. The scarce experimental statistics 
in this energy range does not allow to make the comparison with better accuracy. 



4 Monte-Carlo generator for production of pion pairs 

The same ideas as for muons were applied here to construct the master formula to the processes 
e + e" — > 7r + 7r~(ri7), K + K~(n / -f), KsKjj^ny), assuming that pseudo scalar mesons are point- 
like objects. The enhanced contributions to the cross section, coming from the collinear region, 
are accounted for by means of SF formalism. The one-loop virtual corrections, radiation of 
soft as well as one hard photon are taken into account in the first order of a exactly. The 
effects of the vacuum polarization are not included into the formulae presented below. In this 
case we deal with the so-called dressed cross section, when the dynamics of the pions strong 
interaction is encoded in the form factor properties. But, the Coulomb interaction in the final 
state should be included into RC to eliminate electromagnetic corrections. 

According to the papers [HJ HE] the boosted Born cross section is given by the expression 

daf e -- + *-( Zl , z 2 ) a 2 (Y? - ml/e^ (1 - cjj^g^ 

dfi! '4S z\z\ Zl+Z 2 + ( Z 2- Zl )(l- m 2/( e 2y2))-l/2 Cl ' ^ > 

where z\$, are the energy fractions of the electron and positron after radiation of photon jets 
in the collinear region, \F w (sziZ2)\ 2 is a pion form factor squared, C\ = cos9 1 , Q x is a polar 
angle between momentum of negative pion and electron beam direction. The energy fractions 
Y12 of the final pions and the polar angle of the positive pion, # 2) can be found from the same 
kinematic relations as for muons. 

The charge-even part of the cross section due to radiation of soft and virtual photons [TT^fTR] 
can be written in a convenient way as in jU]: 



dQi dOx 7T 



(A + A n ) 



A, = (L - 1) l„f + \(L - 1) + £ - i, ^, = (i±^,ni±| - l) l„f + AT „„. (25) 
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The expression for the quantity K£ ven can be found in P ITS]. 
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The charge-odd part of the differential cross section is the interference result of the Born 
amplitude with those describing box-type diagrams and soft photons emission by electrons 
and pions |19j . According to jU] the expression for the charge-odd part has the following form: 
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The cross section of pion pair production accompanied by hard photon emission can be pre- 
sented in the following form jHj: 
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The quantity -R^ard * 7 contains the terms describing initial and final state radiation and 
their interference: 
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Mandelstam variables and x±> X± are defined as for electrons. The expression for the master 
formula, describing the process of pion pair production with two compensators, has a similar 
form as for muons and reads 
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(32) 



where = 7r 2 /6 — 1/4 + i^^ en (s, #i) + iQ dd (s, #i), #i is a polar angle of negative pion in the 
center-of-mass system. As well as for muons, the integration limits with energy in the first 
term in Eq 13*21 were again divided in two parts. Two terms describing one photon jet radiation 
are merged with two compensators. Comparison with the BABAYAGA [201 generator was 
performed. The theoretical accuracy of the formulae, used in the BABAYAGA code, is about 
~1%. The current version of the BABAYAGA code (3.5) doesn't include the FSR and so this 
term was removed from our code (just for comparison). The difference of the cross sections 
calculated by the MCGPJ generator and BABAYAGA is shown in Fig[*[7| with the same 
selection criteria as for Bhabha scattering events. A systematic shift between cross sections is 
on average about 1% in agreement with the BABAYAGA code precision, but for the lowest 
and highest energies the agreement becomes worse. 

The distributions of pion, muon and electron pairs as a function of momentum are pre- 
sented in FiglTBI at the c.m.s. energy of 390 MeV for experimental and simulated events. 
Momentum and angle resolutions, decays in flight, interaction with the detector material and 
other important factors were smeared with the parameters of simulated events to create events 
as close as possible to the real ones. The histograms for each type of particles were fitted by 
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two Gaussian functions. Their relative weights and widths were the free parameters of the fit. 
Perfect agreement between experiment and simulation one can see. 

The enveloping curve describes pretty well the shape of histograms both at the peaks and 
at the tails. It permits to determine the number of events inside each histogram and to 
estimate the amount of the muons and electrons under the pion peak and thereby to extract 
the systematic error due to events separation procedure. 

The shape of histograms peaks of the simulated events is not described well, if the MC 
generator, based on the formulae in the first order in a, is used. The shape of the histogram 
peak is mainly driven by the emission spectrum of soft photons and the apparatus resolution. 
Thus, the number of events in the tail area is determined by the peak shape and hence, the 
approach with photon jet radiation is absolutely necessary. 

The MC generator simulating production of charged kaons is created similarly to that for 
pions. The pion mass m n and form factor should be replaced in the above expressions by the 
kaon ones. The cross section being multiplied by the exact Coulomb factor will interpolate 
the energy dependence of the cross section from the threshold production to the relativistic 
region. The exact expression for the Coulomb factor was obtained by Sommerfeld-Sakharov 
and reads 

„, , z , 1i\a 

f(z) = - — - z 2, z = 33 

1 — exp(— z) v 

where v is a relative velocity of koons. The term z/2 is subtracted because it is already in- 
cluded in the 0(a) corrections to final state. 

The MC generator simulating production of neutral kaons is significantly simpler since there 
is no Coulomb interaction and photon emission in the final state. The quantity -R^ard ^ XlXs7 (|3*T]) 
consists of one term which describes the initial state radiation only and the value K is equal 
to vr 2 /6 - 1/4. 



5 Summary and concluding remarks 

The MC generator for the processes e + e~ — > e + e~, 7r + 7r - , K + K~ and KlKs is de- 

scribed in detail. An extended treatment of radiative corrections is implemented in the genera- 
tor to get a high level of theoretical precision. The current version of the program, Monte-Carlo 
Generator Photon Jets (MCGPJ), includes radiation effects in the first order in a exactly. The 
corrections deal with radiation of hard photon decomposed into the three parts which describe 
initial and final state radiation and their interference. All terms in the matrix elements which 
are proportional to the muon or pion mass squared are kept. The enhanced contributions 
coming from the collinear region are accounted for by means of the SF formalism. As a result, 
the theoretical accuracy of the cross sections with RC is estimated to be at ~ 0.2% level. It 
is better by at least a factor of two compared to the accuracy 0.5 — 1% achieved in earlier 
papers. Comparison with the well known codes BHWIDE, KKMC and BABAYAGA shows a 
satisfactory level of agreement for many distributions simulated by the generators. 

The shape of the distributions with acollinearity angles A8 and A0 agrees with the CMD-2 
experimental data. The double ratio of the number of muon events to that of electrons divided 
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by the ratio of the theoretical cross sections was found to be 0.986 ± 0.014. The deviation from 
unity is —1.4 ± 1.4%. This is the first direct comparison of the experimental cross sections 
with the theoretical calculation at the accuracy ~ 1%. The comparison of momenta distribu- 
tions at the lowest energy point shows that only the simulation with radiation of photon jets 
describes the experimental spectra pretty well. Relying upon the above sketched review, the 
main conclusion is that theoretical predictions aiming at a 0(0.1%) precision must include 
contributions of both exact 0(a) terms and all higher order 0(a n L n ) corrections. 

The theoretical uncertainties of the cross sections with RC are determined by the unac- 
counted higher order corrections and they are estimated to be at ~ 0.2% level. Below, the 
main sources of uncertainties in the current formulae are listed: 

• The weak interaction contributions are omitted in our approach. The numerical estima- 
tions show that for energies 2e < 3 GeV these contributions do not exceed 0.1%. 

• A part of the second order next-to-leading radiative corrections proportional to (a /tt) 2 L ~ 
10 -4 were omitted. Among these contributions are: the effect due to double photon emis- 
sion (one inside and one outside of the narrow cone); emission of soft or virtual photon 
simultaneously with radiation of one hard photon at large angles. Even if we assume 
that a coefficient in front of these terms will be of order of ten, their contribution can 
not exceed ~ 0.1%. 

• The next source of uncertainties is related to the calculation of the hadronic vacuum 
polarization contribution to the photon propagator. Numerical estimations show that a 
systematic error of hadronic cross sections of about 1% changes the cross section by less 
than ~ 0.04%. 

• The uncertainty of about 0.1% is related to the theoretical models which are used to 
describe the energy dependence of the hadronic cross sections. 

• In Ref. [21] it was concluded that a combined effect of all the parametrically enhanced 
0(a 2 ) corrections can be numerically limited by 2.0 x 10 -4 for near threshold production. 
The magnitude of this contribution slowly decreases with the final particle velocity (3 
and therefore these corrections are beyond the intended accuracy. 

• The last source of uncertainty is mainly driven by the collinear kinematic approxima- 
tion. Several terms proportional to (a/n)^ and (a/7r)(l/7# ) 2 were omitted. Numerical 
estimations show that the contribution of these factors is about ~ 0.1%. 

Considering the uncertainty sources mentioned above as independent, we can conclude that 
the total systematic error of the cross sections with RC is less than 0.2%. An indirect confir- 
mation of the correct evaluation of the accuracy is the comparison of cross sections with RC 
calculated in the first order of a only. The corresponding difference does not exceed 0.2%. 
From that follows that higher order enhanced contributions, coming from collinear regions 
with emission of two and more photons, contribute to the cross section by amount of ~ 0.2% 
only for our selection criteria. Since the accuracy of this contribution is certainly known better 
than 100%, the systematic theoretical uncertainty for the cross sections with RC is ~ 0.2% 
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Figure 2: The cross section dependence on 
the auxiliary parameter Ae. Parameters 
and selection cuts are given in text. 
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Figure 3: The cross section dependence 
on the auxiliary parameter 9 after inte- 
gration over the remaining kinematic vari- 
ables. 
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Figure 4: The acollinearity polar angle AO 
distribution. The solid line-MCGPJ code, 
the dashed line-BHWIDE. 



Figure 5: The acollinearity azimuthal 
angle A0 distribution. The solid line - 
MCGPJ code, the dashed line - BHWIDE. 
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Figure 6: The distribution of events as a 
function of the missing energy radiated by 
electrons and positrons. The solid line - 
MCGPJ code, the dashed line - BHWIDE. 
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Figure 7: The relative difference of 
cross sections calculated by the MCGPJ 
code and BHWIDE as a function of the 
cm. energy. 
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Figure 8: The relative difference between 
cross sections calculated by the MCGPJ 
code and BHWIDE versus the acollinear- 
ity angle |A0|. 



Figure 9: The relative difference between 
cross sections calculated by the MCGPJ 
code and the generator based on Ref. [7] 
versus the acollinearity angle \A8\. 
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Figure 10: The acollinearity angle A8 dis- 
tribution in the scattering plane. Solid line 
- simulation, histogram - experiment. 
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Figure 11: The acollinearity angle A0 dis- 
tribution in the azimuthal plane. Solid line 
- simulation, histogram - experiment. 
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Figure 12: Two-dimensional plot of the 
simulated events (MCGPJ). The points in 
this plot correspond to the electron and 
positron energies. The influence of the 
condition A8 < 0.25 rad can be seen from 
an ark-like smooth border. 



Figure 13: Two-dimensional plot of the 
simulated events. The generator is based 
on j7j. The points in this plot correspond 
to the electron and positron energies. The 
condition A8 < 0.25 rad divides the plot 
into two parts - with and without events. 
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Figure 14: The difference between cross 
sections as a function of the cut imposed 
on the transverse momenta of final parti- 
cles. 
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Figure 15: The relative difference be- 
tween the cross sections calculated by 
the MCGPJ code and KKMC versus the 
c.m.energy. 



Figure 16: The ratio of the number of se- 
lected muons to that of electrons divided 
by the ratio of the corresponding theoret- 
ical cross sections. 
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Figure 17: The relative difference be- 
tween the cross sections calculated by the 
MCGPJ code and BABAYAGA versus the 
beam energy. 



Figure 18: The distribution of pion, muon 
and electron pairs as a function of the mo- 
mentum. The upper curve represents a 
common fit, bottom curve - mainly cosmic 
ray background. 
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